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There has been much recent interest in directly measuring the electric dipole moments (EDM) 
of the proton and the electron, because of their possible importance in the present day observed 
matter/antimatter imbalance in the universe. Such a measurement will require storing a polarized 
beam of “frozen spin” particles, 15MeV electrons or 230 MeV protons, in an all-electric storage ring. 
Only one such relativistic electric accelerator has ever been built—the 10 MeV “Electron Analogue” 
ring at Brookhaven National Laboratory in 1954; it can also be referred to as the “AGS Analogue” 
ring to make clear it was a prototype for the AGS proton ring under construction at that time at 
BNL. (Its purpose was to investigate nonlinear resonances as well as passage through “transition” 
with the newly-invented alternating gradient proton ring design.) By chance this electron ring, 
long since dismantled and its engineering drawings disappeared, would have been appropriate both 
for measuring the electron EDM and to serve as an inexpensive prototype for the arguably more 
promising, but ten times more expensive, proton EDM measurement. 

Today it is cheaper yet to “resurrect” the Electron Analogue ring by simulating its performance 
computationally. This is one purpose for the present paper. Most existing accelerator simulation 
codes cannot be used for this purpose because they implicitly assume magnetic bending. The 
new UAL/ETEAPOT code, described in detail in an accompanying paper, has been developed for 
modeling storage ring performance, including spin evolution, in electric rings. Illustrating its use, 
comparing its predictions with the old observations, and describing new expectations concerning spin 
evolution and code performance, are other goals of the paper. To set up some of these calculations 
has required a kind of “archeological physics” to reconstitute the detailed Electron Analogue lattice 
design from a 1991 retrospective report by Plotkin[l] as well as unpublished notes of CourantJj 
describing machine studies performed in 1954-55. 

This paper describes the practical application of the ETEAPOT code and provides sample re¬ 
sults, with emphasis on emulating lattice optics in the AGS Analogue ring for comparison with the 
historical machine studies and to predict the electron spin evolution they would have measured if 
they had had polarized electrons and electron polarimetry. 

Of greater present day interest is the performance to be expected for a proton storage ring 
experiment. To exhibit the ETEAPOT code performance and confirm its symplecticity, results are 
also given for 30 million turn proton spin tracking in an all-electric lattice that would be appropriate 
for a present day measurement of the proton EDM. 

The accompanying paper “ETEAPOT: symplectic orbit/spin tracking code for all-electric storage 
rings” [5 documents in detail the theoretical formulation implemented in ETEAPOT, which is a new 
module in the Unified Accelerator Libraries (UAL) environment [5]. 

PACS numbers: 14.20.Dh, 29.20.Ba, 29.20.db, 29.90.+r, 42.25.Ja 


I. INTRODUCTION 

a. Motivation for electric storage ring “traps” 
for electrons or protons. The U.S. particle physics 
community has recently updated its vision of the future 
and strategy for the next decade in a Particle Physics 
Project Prioritization Panel (P5) Report. One of the 
physics goals endorsed by P5 is measuring the EDM of 
fundamental particles (in particular proton, deuteron, 
neutron and electron). 

Since Standard Model EDM predictions are much 
smaller than current experimental sensitivities, detec¬ 
tion of any particle’s non-zero EDM would signal dis¬ 
covery of New Physics. If of sufficient strength, such 
a source could support an explanation for the observed 
matter/antimatter asymmetry of our universe. A pro¬ 


ton EDM collaborationH has proposed a storage ring 
proton EDM measurement at the unprecedented level 
of 10 _29 e-cm, an advance by nearly 5 orders of mag¬ 
nitude beyond the current indirect bound obtained using 
Hg atoms. 

The proposed EDM measurement is based on the ac¬ 
cumulation of tiny “wrong-plane” (i.e. “right” for EDM, 
“wrong” for MDM) spin precessions that can be ac¬ 
counted for only by a non-zero EDM. Polarized beams 
can be frozen (for example in longitudinally-polarized 
states) in storage rings containing appropriate combi¬ 
nations of electric and magnetic bending elements; the 
relative amount depends on the magnetic dipole moment 
of the particles being stored. Only for a few particles, 
of which the electron and the proton are the most im¬ 
portant, can the spins be frozen in purely electric rings. 
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This is highly advantageous, since such rings support 
beams circulating both clockwise and counter-clockwise, 
permitting measurements for which important system¬ 
atic errors cancel. The electric dipole moment (EDM) 
of the proton is known to be so small that, in order to 
reduce systematic errors enough to measure it in a stor¬ 
age ring requires measuring differences between counter- 
circulating beams. This can be done sequentially or, to 
provide differential beam position precision at the cost 
of beam-beam complications, with simultaneous circu¬ 
lating beams. (Direct colliding beam-beam effects are 
calculated to be negligible.) 

A comparably important advantage of electric bend¬ 
ing is that the absence of intentional magnetic fields will 
reduce the presence of unintentional (radial) magnetic 
field components, which are expected to be the domi¬ 
nant source of spurious precession, mimicking the EDM 
effect. 

To freeze the spin procession, conventional storage ring 
bending magnets are replaced with the corresponding 
electric elements. Though the circulating particles are 
constantly being attracted to the inner electrode, their 
centrifugal force causes them to circulate indefinitely in 
a more or less circular orbit. This establishes the storage 
ring as an electric “trap”, albeit for an intense moving 
bunch rather than for a few slow particles. 

In a storage ring EDM measurement, bunches of longi¬ 
tudinally polarized protons will circulate for “long” time 
intervals such as 1000 seconds. Because of inevitable pa¬ 
rameter spreads, individual particle spins will precess dif¬ 
ferently and, after a spin coherence time SCT, the beam 
polarization will have been attenuated (due to decoher¬ 
ence) to a point where the EDM precession rate has be¬ 
come unmeasurably small. The wrong-plane polarization 
difference between early and late times, when ascribed to 
the torque of the bend field acting on the electric dipole 
moment, will provide a measurement of the proton EDM. 

For measuring EDM’s the accumulation and measure¬ 
ment of small effects requires analysis, mitigation and 
control of various systematic errors. Issues to be studied 
with ETEAPOT include: 

Radial B-field: since the torque due to any residual 
radial B-field acting on the magnetic dipole moment 
(MDM) mimics the EDM effect, this is expected to be 
the dominant systematic error. 

Geometric phase: some EDM-mimicking precessions 
would average to zero except for the non-commutativity 
of 3D rotations. 

Non-radial E-field, vertical quad misalignment, RF 
cavity misalignment etc.: these cause the closed orbit 
to deviate from design. The systematic errors they cause 
cancel to “lowest order” but higher order effects need to 
be investigated. 

Polarimetry: realistic beam distributions need to be 
used to identify and reduce left-right asymmetry bias, 
which is another source of systematic error. 

b. Initial storage ring simulation tasks. Initial 

tasks for the ETEAPOT code are to simulate the per¬ 


formance of arbitrary electric storage rings. As with any 
storage ring, these tasks (which now include also spin 
tracking) are 

1. Evaluation of linearized lattice functions such as 
tunes, the Twiss, o, yd, and 7 functions, dispersion 
function, and spin tunes, as well as chromatic de¬ 
pendence of these functions, and their sensitivity 
to imperfections. 

2 . Short term tracking, for confirmation of the lattice 
functions, for determining dynamic aperture, and 
for investigating the performance of analytically- 
derived spin decoherence compensation schemes. 

3. Guaranteed to be stable , long term tracking. The 
EDM experiment requires SCT to be 1000 s or 
longer. During this time every particle executes 
about 10 9 betatron oscillations. The code is re¬ 
quired to exhibit negligible spurious growth or de¬ 
cay for this interval of time. 

Compared to magnetic bending, electric bending compli¬ 
cates these tasks. This requires the ETEAPOT treat¬ 
ment to differ from the well-established TEAPOT treat¬ 
ment. As implemented in UAL/TEAPOT, lattice func¬ 
tion determinations use a truncated power series formal¬ 
ism that has not (as yet) been established for electric 
rings. Translating the arbitrary order formalism from 
magnetic to electric elements is a huge task that has 
only just begun. To cover all three storage ring simu¬ 
lation tasks on the time scale required for EDM planning 
we have, therefore, had to proceed along more than one 
track. 

One track starts by implementing exact (and therefore 
exactly symplectic) tracking in the inverse square law po¬ 
tential electrical field between spherical electrodes. The 
optimal electric field shape will differ from this, however. 
This makes it necessary to introduce artificial virtual 
quadrupoles in the interiors of bend elements, in order 
to model deviation of the actual electric field away from 
its idealized representation. Real quadrupoles are typi¬ 
cally also present in the lattice, for example to alter the 
lattice focusing, or to control dispersion. Thin sextupoles 
(present for example to adjust chromaticities or to com¬ 
pensate spin decoherence effects) and other thin multi¬ 
pole elements are also allowed. Causing only “kicks”, 
these thin elements, either real or virtual, also preserve 
symplecticity. Now complete, this track allows tasks ( 2 ) 
and (3) to be completed for arbitrary lattices. 

Since (linearized) transfer matrices are not used, the 
ETEAPOT code cannot be used to extract Twiss func¬ 
tions directly, as is commonly done in conventional 
Courant-Snyder accelerator formalism. However, tunes 
and Twiss functions can be obtained (to amply satisfac¬ 
tory accuracy) using FFT and MIA (model-independent 
analysis). Now complete this code also provides transfer 
matrix determination which, in turn, allows determina¬ 
tion of all the lattice functions required to complete task 

(i). 
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Implementation of the UAL truncated power series ap¬ 
proach for electric rings has just been begun. Its progress 
is contingent on obtaining government funding[7 U 

After comparing ETEAPOT simulation results to re¬ 
sults measured on the AGS Analogue ring this paper 
briefly describes simulation of spin evolution in a proton 
EDM storage ring. 


II. THE AGS ELECTRON ANALOGUE RING 

Of the more than 100 relativistic accelerators ever 
built, only one has used electric rather than magnetic 
bending. It was the AGS Electron Analogue machine at 
BNL. Curiously it was also the first ring ever to use al¬ 
ternating gradient (AG) focusing. (The Cornell 1.1 GeV 
alternating gradient electron ring was commissioned at 
more or less the same time and the BNL AGS ring itself 
somewhat later.) 

Along with using electrons instead of protons , and lim¬ 
ited by achievably high electric field, cost minimization of 
the AGS Analogue led to the choice of 10 MeV maximum 
electron energy, 4.7 m bending radius, 6.8 MHz RF fre¬ 
quency, and 600 V RF voltage. These optimization con¬ 
siderations are very much the same as will be used to fix 
the parameters of a frozen spin proton ring (which is ten¬ 
tatively expected to have a bend radius of about 50 m) [ 61 . 
The UAL/ETEAPOT code (documented in the accom¬ 
panying paper) was developed with this application in 
mind. The present paper gives initial results. 

Starting from fragmentary documentation, the pa¬ 
per starts by reverse engineering the AGS Ana¬ 
logue lattice and producing a lattice description file 
E_AGS_Analogue . sxf, in the format needed for process¬ 
ing in the Unified Accelerator Libraries (UAL) environ¬ 
ment. Results obtained using ETEAPOT are then com¬ 
pared with measurements performed on the ring at BNL 
in 1954-55. 

By chance, the magic kinetic energy for freezing elec¬ 
trons, which is 15Mev, is not very different from the 
10 MeV of the AGS Analogue ring. So that ring could 
have been used to measure the electron EDM just by in¬ 
creasing the electric bend field by a factor of 1.5. Morse 0 
and others have suggested building such a ring for this 
purpose. The ETEAPOT code can therefore be used to 
to simulate an electron EDM measurement using a ring 
whose successful performance as a storage ring is all but 
guaranteed by the successful operation of an “identical” 
ring in 1955. 15 MeV is a very convenient electron energy 
and high quality electron source would be available, for 
example as described by Bazarov[8]. 

The present paper describes our “resurrection” of AGS 
Analogue ring from historic BNL documentation, and 
simulates its performance with codes intended for the 
proton EDM experiment. 

A quite superficial (day long) search of the BNL li¬ 
brary and the Accelerator-Collider report library found 
one quite extensive report, produced retrospectively in 


1991 by Martin Plotkinfl], and a few ancient reports de¬ 
scribing machine studies results. Also a report privately 
communicated from Ernest Courant [2| contains experi¬ 
mental data to be simulated in this paper. 

Before starting on this project, one of the authors, RT, 
benefited from three brief but valuable meetings with 
Ernest Courant, perhaps the father, or at least one of 
the parents, of the Electron Analogue ring, inquiring 
about his recollections concerning the dynamics and per¬ 
formance of the ring. This contributed to our reconstruc¬ 
tion. Bill Morse [3] reports having had a similar conversa¬ 
tion with Courant a few years earlier. Bill recalls asking 
Ernest whether, back in 1953, he (Courant) understood 
the difference between electric and magnetic focusing. 
Ernest’s replied “of course” in his always kindly, but in 
this case somewhat exasperated, tone of voice. 

c. Historical BNL documents. The “Concep¬ 
tual Design Report” for the AGS Analogue electron ring 
was a four page letter, dated August 21, 1953, from BNL 
Director Haworth to the A.E.C. (predecessor of D.O.E.) 
Director of Research Johnson, applying for funding. The 
letter is reproduced in its entirety in Figure [7] in the ap¬ 
pendix. As brief as it is, this letter along with hints from 
Plot kin and Courant, includes everything needed to re¬ 
construct the ring, as shown in Figure [l] By 1955 the 
ring had been approved, built, and commissioned, and 
had achieved its intended purpose. 

Following the Haworth letter in the appendix are two 
other especially informative figures from the paper by 
Plotkin. Figure [8] is especially useful for visualizing the 
physical layout of the AGS Analogue ring and its vacuum 
system. Figure [9] is a more polished version of Figure [3] 
which plots tune scans actually performed on the AGS- 
Analogue ring and reported by Ernest Courant in a July 
28, 1955 BNL technical report [2]. 

In these plots, points of observed beam loss and ob¬ 
served beam disruption in the ( Q x ->Qy ) tune plane are 
correlated with expected resonances. Beam loss occurs 
on integer resonances, beam disruption occurs on half in¬ 
teger resonances. The Courant report on, and analysis 
of, data collected in machine studies less than two years 
after the submission of the ring funding proposal, would 
certainly deserve an A+ grade by modern machine stud¬ 
ies grading standards. 

The axes are voltages (proportional to focusing 
strengths) applied to the tune-adjusting quadrupole fam¬ 
ilies. Short heavy lines indicate regions with no beam sur¬ 
vival (presumably due to integer resonance). Dots indi¬ 
cate points reported by Courant as “showing the charac¬ 
teristic double envelope of the oscilloscope pattern, some¬ 
times accompanied by beam loss”. (These were presum¬ 
ably due to “very narrow” 1/2 integer resonance). The 
nominal central tunes values are ( Q x ,Qy ) = (6.5, 6.5). 
Stop bands due to the eightfold lattice symmetry are also 
shown. 
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FIG. 1. The 1955 AGS-Analogue lattice as reverse engineered from available documentation—mainly the 1953 proposal 
letter from BNL Director Haworth to AEC Director of Research Johnson. Individual “lenses” (as they were referred to 
in the original documentation; in modern terminology they would be referred to as “combined function bends”) are shown 
on the left; vertical lines indicate locations of thin quadrupoles present in the lattice description to represent the focusing 
effect of the electric lens elements; Elements are identified by ad hoc labels assigned in the reconstruction. They are not 
obtained from the original lattice design, but are needed for modern day lattice description files. Lattice description files are 
avaiable in various formats: E_AGS_Analogue2.xml gives the pristine design with all parameters in algebraic (or function of 
algebraic) form. E_AGS_Analogue2. adxf is the same, but all parameters values, though ideal, converted to numerical values. 
E_AGS_Analogue2. sxf contains the fully-instantiated lattice description (with parameters from otherwise identical elements 
allowed to be individualized). Most UAL lattice files are available in all these forms. 


III. CURRENT DAY SIMULATION OF 1955 
MACHINE STUDIES TUNE PLANE SCAN 

Our AGS Analogue lattice reconstruction is shown in 
Figure [l] The Courant tune plane plot is shown in Fig¬ 
ure [2] It is to be compared with a similar plot, simulated 
by TEAPOT and shown in Figure [3J As far as we know 
this code and this lattice representation are completely 
equivalent Courant’s model and analysis in 1955. 

In the TEAPOT tune plane plot, boxes indicate points 
on integer resonance boundary curves of the stable di¬ 
amond centered on nominal tune values (Q x ,Qy) = 
(6.5,6.5). Points lying on 1/2 integer resonance lines 
are indicated by dots. Superperiodicity (eightfold pe¬ 
riodicity) causes resonances with Q x = 8 or Q y = 8 indi¬ 
cated by broad blank bends bounded by narrowly-spaced 
lines. Courant refers to these as “stop bands”. Hor¬ 
izontal/vertical axes (Q + ,Q _ ) are “electrode voltages 
on the quadrupoles in odd/even-numbered tanks”. (For 
the meaning of “tank” see Figure [8j) The quadrupole 
strength coefficients for variable quadrupoles (Q + ,Q - ) 


were determined empirically to match the central tunes. 
This means that absolute focusing strength scales are not 
checked. Otherwise there are no significantly adjustable 
empirical lattice parameters. 

The AGS Analogue ring provides only a coarse test 
of ETEAPOT since, for strong focusing, the change 
from electic and magnetic bending is quite minor. 
This is illustrated in Figure [4] which shows that the 
(Q x , Qy)=(6.5,6.5) tunes in the AGS Analogue ring are 
high enough that the tune plane structure is quite insen¬ 
sitive to whether the bends are treated as magnetic or 
electric. 

Treating the difference perturbatively, the effects of 
changing from magnetic to electric bends are inversely 
proportional to /?’s, which are both 6.5 in this case. For 
the eventual proton EDM ring the vertical tune has to be 
reduced from Q y = 6.5 by at least a factor of ten which 
completely invalidates any such perturbative estimation 
and requires radical adjustment of quadrupole strengths. 
Reducing Q y from 6.5 to 2.25 has been straightforward 
but, as the electric focusing became increasingly impor- 
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tant, to decrease Q y further will require substantial de¬ 
sign effort. 

Spin evolution in the AGS Analogue ring is shown in 
Figure [5] and described in the caption. From these results 
we are confident in out understanding of electric rings 
and of our ability to simulate their performance using 
ETEAPOT. 
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FIG. 2. Tune plane resonance diagram measured during machine studies at the AGS Analogue ring and reported by Courant[2|. 
A more polished and complete version of this plot, copied from PlotkinpQ is shown in Figure [9] but with minor re-labelings. 



Q' = 40000 qklx-7811 


FIG. 3. Tune plane resonance diagram as calculated by TEAPOT, with bend elements treated as magnetic. There is good 
qualitative and quantitative agreement with the Courant data shown in Figure [2] 
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Beta_x functions AGS_Analogue electric m=-l Beta_y functions AGS_Analogue 




position [m] position [m] 


FIG. 4. Comparison for AGS_Analogue2. sxf lattice with bends treated as electric m = 1 or m = -1 or magnetic. f3 x is plotted 
on the left, j3 y on the right. As explained in the text, because of the strong focusing in the AGS Analogue lattice, switching 
from magnetic to electric bending causes only minor changes to the lattice functions. 
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FIG. 5. The upper Hots show the s x and s y spin components of a single electron, with initial horizontal amplitude xtyp=le6m, 
in the AGS Analogue ring. (The 6D initial conditions are given above the plots.) The horizontal axis “split bend index” 
increases by 1 at each bend, of which there are 640, so about 59 turns are shown. The lower right figure shows longitudinal 
phase space evolution for three off-energy particles, with energy offset de plotted against time offset —ct. One sees that almost 
exactly one revolution is completed during this time. The synchrotron tune is therefore approximately 1/59. The lower left 
plot plot shows the de energy deviation during one synchrotron oscillation period. 




















IV. LONG TERM TRACKING IN A PROTON 
EDM STORAGE RING TRAP 

Plots so far have simulated the evolution of lOMeV 
electrons in the AGS Analogue electron ring. These 
roughly indicate the performance to be expected in a 
modern day electron EDM storage ring measurement. 
Of more immediate interest is the expected behavior of 
230 MeV protons in a significantly larger, but still all¬ 
electric proton EDM storage ring. 

A single particle tracking example, taken from a re¬ 
cent proton EDM study, is illustrated in Figure [6j The 
particle orbit and spin components are tracked around a 
prototype proton EDM storage ring for 33 million turns 
using ETEAPOT. As stated earlier the tracking is exact 
and no artificial “symplectification” is applied. Any spu¬ 
rious damping or anti-damping of the spatial orbit over 
the full run is less than ten percent over the 33 million 
turns. The (unit-magnitude) spin vector S is initially 
purely tangential (forward) and the vector magnitude 
never changes from its initial value of 1. 

It is characteristic of spin evolution for the transverse 
spin components to change over narrow ranges in syn¬ 
chronism with betatron oscillation and larger but still 
small, ranges in synchronism with synchrotron oscilla¬ 
tion. The upper left graphs of Figure [6] show this, along 
with a sinusoidal fit with the parameters shown. These 
precessions are not expected to contribute significantly to 
spin decoherence. But any systematic growth over mil¬ 
lions of turns will eventually lead to beam decoherence 
and limited spin coherence time. 

An important computational task in planning to mea¬ 
sure the proton EDM is to determine the spin coherence 
time SCT of the circulating beam. This is the time after 
which inevitable spreads in beam parameters will have 
attenuated the beam polarization significantly due to de¬ 
coherence in the spin propagation. As long as both S x 
and S y magnitudes remain small, S remains in the for¬ 
ward hemisphere, and decoherence is suppressed. 

One sees from the graphs of Figure [6] that both S x 
and S y gradually deviate from zero. But, for the partic¬ 
ular particle being tracked, these deviations are limited 
to small values. In particular, since S% + Sy remains 
always much less than 1 in magnitude, S z therefore re¬ 
mains not much less than 1. So the beam polarization re¬ 
mains always in the forward hemisphere. Depolarization 
is therefore unimportant for particles of amplitude com¬ 
parable to, or smaller than, this particle’s. The tracking 
is therefore consistent with the SCT value exceeding the 
33 million turns shown. 

The lattice investigated in Figure [6] is a “Mobius lat¬ 
tice” in which the horizontal and vertical betatron oscilla¬ 
tions interchange on every turn. This strongly suppresses 
spin depolarization because precession in horizontal and 
vertical oscillation phases tend to cancel on a turn-by¬ 
turn basis. By fine tuning this cancellation it is antici¬ 
pated that extremely long SCT values can be obtained. 

This will leave beam energy spread as the dominant 


source of beam depolarization. Fringe fields at bend el¬ 
ements are one such significant source of spin precession 
slewing. In this respect spin tracking is more sensitive 
than orbit tracking. Fringe fields are not very important 
in orbit calculations. Treating bends as “hard edged” 
mainly causes small tune shifts which are not very im¬ 
portant, since both tunes are always set operationally 
using spectrum analysis. But, because of the spin pre¬ 
cession sensitivity, in ETEAPOT, fringe fields are treated 
as linear ramps of length comparable with the gap be¬ 
tween the electrodes, and not by hard edges. In any case 
synchrotron oscillation averaging also strongly suppresses 
this fringe field depolarization. This places further de¬ 
mands on longitudinal particle tracking, which is always 
delicates, because nonlinearity of synchrotron oscillations 
contributes to spin decoherence. 

The particle revolution period is about one microsec¬ 
ond so the plots shown correspond to about thirty sec¬ 
onds of real clock time in the laboratory. In a laptop 
computer this computation takes a few hours per par¬ 
ticle. The ratio of computation time to real laboratory 
time for a single particle is in the range from one to ten 
thousand. The present ETEAPOT single particle track¬ 
ing approach is sufficient for early design tasks. 

EDM storage ring beams will contain perhaps 10 8 par¬ 
ticles. Adequately precise simulation of the EDM ex¬ 
periment will require fewer than this, but at least thou¬ 
sands of particles to be tracked for significantly longer 
times than shown. This will require heavy paralleliza¬ 
tion of the tracking. With little particle-to-particle in¬ 
teraction, the code is easily parallelizable. For more ad¬ 
vanced tasks, such as investigating polarimeter biases or 
emittance growth caused intrabeam scattering, an effi¬ 
cient and scalable map-based computational approach is 
under consideration. 
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FIG. 6. Evolution of spin components for 33 million turns. Because S z — 1 initially, both S x and S y vanish at x=0, where x is 
turn number. The upper left figure shows S x for the first 1000 turns, as fit by a sinusoid Asin(wx) where x is turn number. 
The upper right figure shows S x for one million turns starting at x = 29 million. The lower left and upper right figures show 
the evolution of S x and S y for the full 33 million turn run. Though these motions are oscillatory the oscillation amplitudes are 
less than the line width. The lattice file is E_pEDM-rtrl-Mobius .RF. sxf. 
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BNL Director Haworth Electric Analogue Proposal Letter to AEC 

August 21, 1953 


Dr. T.H. Johnson, Director 
Division of Research 
U.S. Atomic Energy Commission 
Washington 25, D.C. 

Dear Tom: 

This letter concerns certain aspects of our accelerator development 
program, particularly the proposed electron model. 

As you know, the general development of a very high energy alternating 
gradient synchrotron is proceeding actively at Brookhaven, utilizing operating 
funds allocated to Basic Physics Research. As I explained in my letter of 
August 12, however, these funds are insufficient to carry forward the 
development as rapidly as desirable. Also, there are certain steps which 
should be taken for which the expenditure of operating funds is not 
appropriate. The first and most important of these is the construction of an 
electron model intended to provide final assurance of the technical 
feasibility of the chosen machine and, more importantly, to provide 
information enabling us to design in the most effective and economical manner. 
(We have no doubt of the general feasibility of accelerators of this type.) 

We have given considerable thought to the requirements for such a model 
and to the philosophy which should guide us in designing and building it. In 
the alternating gradient synchrotron, two problems require especially careful 
exploration by extensive calculation and experimental modelling. These are 
the close-spaced resonances in the betatron oscillations and the shift of 
phase stability at intermediate energies. It seems best to study these 
problems with an electron accelerator which would be essentially an analogue 
rather than an exact model. This device should, in our opinion, be designed 
to yield the maximum of orbital data with a minimum of engineering 
complications, especially those not applicable to a final machine. After 
considerable thought we have arrived at a tentative description and list of 
parameters which follow. 

The device would consist of an accelerator having an orbital radius of 
15 feet and an overall diameter including the straight sections, of 
approximately 45 feet; the guide and focussing fields would be electrostatic, 
with electrode shapes as indicated in the sketch (full scale). 


FIG. 7. This figure, continued on the following three pages, is a copy of the original funding request letter from BNL for the 
AGS Analogue ring. 
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Dr. T.H. Johnson Page 2. 


^ Center 
^ ~ - 


Electrons of about 1 MeV energy would be injected from a small horizontal 
Van de Graaff generator (of the 2 MeV type manufactured by the High Voltage 
Engineering Corporation) so that 5% to 6 % frequency modulation would be 
required. 

Use of a reasonably large radius helps the radio frequency and observing 
equipment in frequency range where good techniques exist, and permits high 
n-values which are necessary for strong alternating*gradient focussing, 

(This, and phase transition, will not be modeled in the Cornell machine.) A 
moderate rise rate, consistent with attainable vacuum requirements, still 
permits the use of small, air-cooled amplifier tubes and a heavily loaded low- 
Q rf cavity. 



A tentative list of parameters is: 


Radius of curvature 15 ft 

Over-all diameter 45 ft 

n 200 

No. of periods 37 

No. of straight sections 74 

No. of lenses per period 4 

Length of lens 7.6 in. 

Length of straight section 7.6 in. 
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Field strength (magnetic type) 



at injection 

10.5 

gauss 

at 10 MeV 

74 

gauss 

Field strength (electrostatic type) 



at Injection 

3 

kV/cm 

at 10 MeV 

22 

kV/cm 

Rise time 

.01 

sec 

Phase transition energy 

2.8 

MeV 

Frequency (final) 

7 

me 

Frequency change 

5k 

% 

Volts/turn 

150 

V 

RF power 

about 1 

kw 

No. of betatron wavelengths 

about 6.2 


aperture 

1X1 

in. 

Betatron amplitude for 10' 3 rad. error 

0.07 

in. 

Maximum stable amplitude, synchrotron 

osc,-0,16 

in. 

Radial spacing of betatron resonances 

about 0.4 

in. 

Vac 1 ium requirement about 10" 6 

mm Hg 


Total pow sr requirements will be small and available with existing 
installations, The test shack seems to be a suitable location since the ring 
will be erected inside a thin magnetic shield which can be thermally insulated 
and heated economically. 


We estimate the cost 
in the following table: 

to be approximately 

§600,000, 

distributed as shown 

Model 

Direct 

Overhead 

Total 

Staff S. fit W. 

§135,000 

§ 65,000 

§200,000 

Van de Graaff 

70,000 

- 

70,000 

Other E. & S, 

130,000 

- 

130,000 

Shops 

135.000 

65.000 

200.000 

§470,000 

$130,000 

$600,000 
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It is seen that total direct expense would be approximately $470,000 and 
the overhead assigned by our methods of accounting would be approximately 
$130,000, Since the total Laboratory overhead cost would not actually be 
appreciably affected by the project, the assignment of this overhead would 
result in relieving Laboratory operations of cost in approximately the same 
amount. A small but finite savings to operations would also be affected by 
way of shop charges since approximately one-third of the shop costs represent 
fixed expenses for foremen, tool crib attendants and other helpers. 

The above facts would, of course, reduce the budgetary difficulties in 
which we find ourselves as a result of the present operational ceilings. 
Furthermore, the construction of the model would result is some transfers of 
scientific, technical and shop personnel from the operating program, resulting 
in further savings. We would naturally like to minimize this by adding new 
personnel either in operations or for the construction of the model. Allowing 
for a reasonable compromise in this connection, I believe that our Physical 
Research Program could be carried forward without too serious sacrifices by 
the provision of the $600,000 required for the electron model, the additional 
$100,000 promised for Applied Chemistry and $350,000 additional operating 
funds for Applied Physics. Were all of this forthcoming it would not, in my 
opinion, be necessary to reduce the power level of the reactor but we would 
endeavor to make all possible economies in equipment and supplies and in 
overhead costs and would hold physics personnel to present levels except for a 
few additions in the accelerator development work. 


Sincerely yours, 


UH:DS 

cc: E.L. VanHorn 


Leiand J. Haworth 
Director 




FIG. 8. Electric “leneses” were contained, 10 to a tank, to the 16 vacuum “tanks” making up the full AGS Analogue ring. 
Each tank therefore contained 2| lattice cells, as well as miscellaneous, special purpose quadrupoles, chromaticity-correcting 
sextupoles and octupoles. 
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RESONANCE DIAGRAM vs. QUADRUPOLE VOLTAGES 

DATA TAKEN FROM ELECTRON ANALOGUE JtUlLY 28, 1955 
tWlTH OPTIMUM SETTINGS OP KICK EHS ft 3EXTUPOLE3 



FIG. 9. This figure, from Plotkinjl], with little accompanying explanation, clearly represents the refined analysis of AGS 
Analogue machine studies data like that shown in Figure [2] 
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